Astronomy <fc Astrophysics manuscript no. 7542 


© ESO 2008 


February 1, 2008 





An approximate theory for substructure propagation in clusters 



D. A. Prokhorov 1 ' 2 and F. Durret 1 



o 
o 

< 



43 

6 



1 Institut d'Astrophysique de Paris, CNRS, UMR 7095, Universite Pierre et Marie Curie, 98bis Bd Arago, F-75014 
Paris, France 

2 Moscow Institute of Physics and Technology, Institutskii lane, 141700 Moscow Region, Dolgoprudnii, Russia 
Accepted . Received ; Draft printed: February 1, 2008 

ABSTRACT 

Aims. The existence of dark matter can be proved in an astrophysical context by the discovery of a system in which the 
observed baryons and the inferred dark matter are spatially segregated, such as the bullet cluster (1E0657-558). The 
full descriptions of the dark matter halo and X-ray gas substructure motions are necessary to forecast the location of 
the dark halo from X-ray maps, which can be confirmed by the detection of a galaxy concentration or by gravitational 
lensing. 

Methods. We present an analytical hydrodynamic model to determine the distance between the X-ray and dark-matter 
components and the Mach number of the merger shock. 

Results. An approximate solution is given for the problem of the substructure propagation in merging clusters. A new 
method to predict the position of a dark matter halo in clusters, where there is a separation between the X-ray gas and 
the dark halo, is proposed and applied to the clusters 1E0657-558 and Abell 1763. 
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1. Introduction 

The "unseen" component of the universe (which dominates 
its mass), is a long-standing issue in modern cosmology 
since the need for dark matter was originally pointed out 
by Zwicky (1933). In a typical cluster of galaxies only 2% 
of the total mass is in the form of galaxies, 10-15% is in 
the form of hot X-ray emitting gas and the remainder is 
dark matter. Modified Newtonian dynamics (MOND) is an 
alternative theory that can be used to model different ob- 
servations without assuming the existence of dark matter 
(Milgrom 1983). However, MOND cannot account for all 
the properties of clusters, where the presence of dark mat- 
ter remains unavoidable (Gerbal et al. 1992, Pointecouteau 
& Silk 2006). 

The existence of dark matter can be confirmed by the 
discovery of a system in which the observed baryons and the 
inferred dark matter are spatially segregated. One type of 
object where this separation occurs in is merging galaxy 
clusters. During the collision of two clusters, the galax- 
ies are effectively collisionless particles, while the plasma 
clouds are highly collisional and therefore are slowed by 
ram pressure. If dark matter particles are also collisionless, 
as is widely assumed, any dark matter present in the sys- 
tem would be located near the galaxies. One such merging 
cluster system is 1E0657-558 (z=0.296), the so-called bul- 
let cluster. The cluster 1E0657-558 has two primary galaxy 
concentrations separated by 0.72 Mpc on the sky, a less 
massive (T=6 keV) western subcluster and a more mas- 
sive (T=14 keV) eastern main cluster (Markevitch et al. 
2002). Both concentrations have associated X-ray emitting 
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plasma offset from the galaxies toward the center of the 
system (Clowe et al. 2006). 

Shock waves driven in the intergalactic medium dur- 
ing the merging of galaxy clusters have been observed in 
X-ray imaging and spectroscopy. Chandra observations of 
1E0657-558 have revealed a bow shock propagating in front 
of a bullet-like gas cloud moving away from the core of 
the main cluster. Based on the gas density jump across 
the shock front, Markevitch et al. (2002) derived a Mach 
number of the shock ~3. The inferred shock velocity of 
~4700 km s _1 has been commonly interpreted as the veloc- 
ity of the "bullet" subcluster itself (Markevitch et al. 2002; 
Hayashi & White 2006; Markevitch 2006, among others). 

Numerical simulations have very recently been run by 
several teams and succeed quite well in reproducing the 
main properties of the bullet cluster, both qualitatively 
and quantitatively. For example, Milosavljevic et al. (2007) 
conclude that the halo collision velocity need not be the 
same as the intergalactic gas shock velocity; their simula- 
tion finds that the velocity of the CDM halo is ~ 16% lower 
than that of the shocks. Springel & Farrar (2007) derived 
a shock speed of about 4500 km s _1 but a subcluster ve- 
locity of only 2600 km s _1 in the rest frame of the system, 
so the shock wave propagates faster than the dark matter 
clump. Randall et al. (2007) combine numerical simulations 
with results derived from X-ray, strong and weak lensing, 
and optical observations to place an upper limit on the 
self-interaction cross-section of dark matter per unit mass 
a/m, which rules out most of the range previously invoked 
to explain inconsistencies between the standard collision- 
less cold dark matter model and observations. Bradac et 
al. (2006) have reconstructed the mass distribution of this 
cluster both from weak and strong lensing data, based on 
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multi band high resolution HST ACS images. They con- 
firmed that the total mass in this cluster does not trace 
the baryonic mass, as already found by Clowe et al. (2004, 
2006). 

For the description of the X-ray substructure propaga- 
tion, we are going to use the "piston in a tube" model, 
described in Sect. 2. In accordance with this model, an ad- 
ditional retarding force acts on the X-ray substructure. The 
analysis of substructure propagation is considered in Sect. 
3 for various cluster parameters. 

Several estimates from the X-ray observations indicate that 
around 30% of all rich clusters exhibit substructure on a 
scale of 1 Mpc (Forman & Jones 1994). The number of 
rich clusters with substructure is a measure of the fraction 
of galaxy clusters that have recently accreted a significant 
fraction of their mass (Richstone et al. 1992), which can be 
used to determine flM- 

The dynamical time of the substructure is derived from the 
approach of Gunn & Gott (1972) in Sect. 4. Unless eval- 
uations of the incidence of X-ray substructures are over- 
estimated or substructure lasts significantly longer than a 
gravitational free-fall time in clusters, the value of the frac- 
tional density FIm is larger than 0.5. The retarding force 
has an effect on the X-ray substructure motion and the 
duration of the merger can be longer than a gravitational 
free-fall time because of this force. 

We also apply the "piston in a tube" model to the clus- 
ter Abell 1763, where the X-ray gas substructure does not 
seem to have gone beyond the main cluster center. A new 
method to predict the position of the dark matter halo in 
clusters, where the observed baryons and the inferred dark 
matter are spatially segregated, is proposed. 



2. The "piston in a tube" model 

The discovery of a galaxy cluster in which the observed 
X-ray plasma and the inferred dark matter are spatially 
segregated requires the existence of a retarding force act- 
ing on the X-ray substructure. The observed shock front 
(in 1E0657-558) moving ahead of the X-ray clump shows 
supersonic motion of the gas. 

We develop a simple 1-D "piston in a tube" model to 
derive the analytical expression for the retarding force. The 
retarding force acting on the X-ray clump (piston) in this 
model arises because the pressure of the compressed gas 
ahead of the piston is higher than the pressure of the rar- 
efied gas behind the piston. Different regions of the flow are 
essential for this model: 1) the region of calm gas ahead of 
the piston, 2) the region of shocked gas, 3) the region of 
rarefied gas, and 4) the region of calm gas behind the rar- 
efied gas. The piston is situated at the boundary of regions 
2, 3 and moves towards region 1. The shock is situated at 
the boundary of regions 1,2. 

In regions of calm gas, 1 and 4, the velocity of the gas 
equals 0, the pressure is that of the ambient medium. In 
region 2 (post-shock region), the velocity of the gas equals 
the piston velocity U, the pressure is higher than the initial 
one. In region 3, the rarefaction wave arises; the velocity 
of the gas depends on the distance from the piston and 
the value of the pressure is smaller than the initial one. 
Qualitative sketches of the gas velocity and pressure profiles 
in different regions are plotted in Figs. 1 and 2, respectively. 



Velocity 




Fig. 1. Gas velocity in different regions of the flow. The 
piston moves from left to right. 

Pressure 
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Fig. 2. Gas pressure in different regions of the flow. 

2.1. The pressure in front of the substructure 
(region 2) 

Following Landau & Lifshitz (1959) we calculate the pres- 
sure in front of the substructure (region 2). The shock wave 
arises in front of the substructure (which acts as a pis- 
ton) because the substructure has a velocity higher than 
the sound velocity in the main cluster. At first the posi- 
tion of the shock coincides with the position of the piston. 
Later the shock wave outdistances the piston and a region 
of shocked gas appears between the shock and the substruc- 
ture. The gas pressure in front of the shock wave (region 1) 
equals the initial pressure p\ and its velocity equals zero. 
The gas in the region between the shock and the substruc- 
ture (region 2) has the substructure velocity. The differ- 
ence between the velocities of regions 1 and 2 equals the 
substructure velocity U. 

The conditions on the shock (continuity of matter flux, 
momentum flux, energy flux) can be written as: 

PWi = p 2 v 2 = j (1) 
Pi + P\v\ =P2 + P2V2 (2) 

2 2 

V V 

Y + w i = y + W2 (3) 

where v is the gas velocity in the reference frame of the 
shock, p is the gas density, p is the pressure, and w is the 
specific enthalpy. 

From the first equation we find the difference of the veloc- 
ities: 

vi-V2 = — - — =j-(V 1 -V 2 ) (4) 
Pi Pi 



D.A. Prokhorov & F. Durret: An approximate theory for subcluster propagation 3 

where V is the inverse of the gas density. 2.2. The pressure behind the substructure (region 3) 

Therefore we can express the matter flux . / 

Ine rarefaction wave is formed behind the substructure (re- 

_ v ^ gion 3). The self-similarity solution of the flow is described 

j = — — (5) in Landau & Lifshitz (1959), this solution depends on the 

1 2 single parameter £ = x/t, where a; is a distance, t is a time 



Using the expression vi. 2 = j • Vi,2 and Eq. (2), we can find 



parameter. The equation of entropy conservation is 



ds ds 

Pi+.fV 1 = P2 +.fV 2 (6) d~t +V dx~ = ° (16) 

, 2 where s is the specific entropy, v is the flow velocity. 

One can express j : rp^g derivatives can be expressed as 



•2 P2 - Pi ,^ d £ d 



>' = ^ (7) s-i« (17) 

did 

We can use the Eqs. (5) and (7) together to calculate the d~x ~ ld£ 
velocity U of the substructure: 

Using Eqs. (16) , (17), (18) one can find 
(w-04=0 (19) 
Using the expression vi. 2 = j ■ Vi i2 and Eq. (3), we can find w here s'^ = ds/d^. 



U = V!-V 2 = V(P2-Pi){Vi - V 2 ) (8) 



The Euler equation and the continuity of matter flux equa- 

(20) 
(21) 



tion are 




dv dv 
dt dx 


1 dp 


p dx 


dp dv 
dt 13 dx 


+ v dp - 
dx 



i 2 V 2 i 2 V 2 
Wl + = W 2 + (9) 

Using Eq. (7), we can find 

w l -w 2 = \{V l + V 2 ){p 1 -p 2 ) (10) 

1 We can rewrite these equations: 

The specific enthalpy of polytropic gas is ( v — Qp' + py' = (22) 

Wl 2 = ™* V ™ (11) (v-ZH = - C ^ (23) 

(7-1) 

Consequently from Eqs. (19) and (22), we have = and 
where 7 = c p /c v . s = cons t an d the self-similarity solution is isoentropic. 

After the simple transformation of Eqs. (10) and (11), one Eliminating p' v' £ from these two equations, one can find 
can find the equation for the percussive adiabat thfi condition for thc non _t riv ial solution (Landau & Lifshitz 

1959): 

V2 = (7 + j>i + (7 - 1)P2 n2] 

Vi ( 7 -l)pi + (7 + l)f>2 1 ' Z- V = c ( 24 ) 

o u t-t. t/ * t? /io\ • /o\ • r( Using the Eqs. (22) and (24), we have 

Substituting U 2 from Eq. (12) m Eq. (8) one can simplify 

the expression for the substructure velocity cp^ — pv^ (25) 

which can be integrated to give v: 

V =[<P^P (26) 
J P 



u = iP2 - pi) \j(j-i)p 1+ \j + i)P2 (13) .._ r«p)*> 

For a polytropic gas p = Apt, therefore the sound velocity: For a diabatic processes 



2 dp . „_! TA pT 1 ^ 1 -^ = const (27) 

c 2 = = 7^p 7 1 = 7pU (14) ' v ; 

" The sound velocity is proportional to the square root of the 

One can find thc expression for the pressure in the region temperature, t ere ore 
2 from Eqs. (13) and (14): , x 2/(7-1) 

P = Pol±) (28) 



P2 , 7( 7 +l)[/ 2 iU L ( 7 + l) 2 £/ 2 

— = 1 + £3 + — y 1 + f^2 ( 15 ) Using Eqs. (26) and (28) we have 

where ci is the sound velocity in the first region. 



v = — ^— [dc= (c-ci) (29) 
7-17 7-1 
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where c\ is the sound velocity in the region of the calm gas. 
7-1, , 



c = C\ 

Using Eq. (28), we obtain the density in the region 3: 

2/(7-1) 



(30) 



P3 = Pi 1 - 



7- 



1 v 



2 a 

and the pressure ps = 



(31) 



P3 = Pi 



1 - 



7' 



! 27/(7- 



Cl 



2.3. T/je retarding force 

From Eqs. (15) and (32) the retarding force F ■■ 
is: 



(32) 



{P3~P2)S 



3. Simulation of substructure propagation 

3.1. Merger initial velocity 

During the merger of two clusters of galaxies, the bary- 
onic component, feeling the gravitational potential cre- 
ated by the dark matter components, moves supersonically 
(Gabici & Blasi 2003). Let us assume that two clusters are 
completely virialized structures with masses Mi and M 2 
(M 2 < Mi, e.g., the cluster 1E0657-558). To approximate 
the merger velocity of the system, we use the assumption 
(Bcrrington & Dermer 2003) that the cluster of mass M 2 is 
a point mass. The main cluster accretes another cluster that 
falls from a turnaround radius Rt a of the system, where the 
two subclusters are supposed to have zero relative velocity. 
The exact value of the turnaround radius does not strongly 
affect the collision velocity as long as it is large and the 
infall velocity approaches free-fall from infinity. Using the 
energy conservation law, one can find the velocity of the 
subcluster with mass M 2 at the virial radius R V i r i of the 
massive cluster, which has the mass M±: 




GMi GM\ 



Rt 



v = 



2GM\ 2GM 1 



R,, 



Rta 



The definition of the virial radius is 



Let 7 = 5/3 and u = U/3ci, we can rewrite Eq. (33) using 
these definitions 



3Mi 



MPc 



1/3 



(36) 
(37) 

(38) 



F(u) = piS ((1 - u) 5 - 1 - lOu - 5m vT" 



Au 2 



If u <C 1, one can calculate the Taylor series 
G{u) = -20 Pl Su 



(34) 



(35) 



where A c is the density contrast for the formation of the 
cluster, S7jvf is the matter density fraction and the critical 
density is 



Pc 



3H 2 
8^G 



(39) 



The functions F/faS) and G/{p x S) are shown in Fig. 3. 3.2. Solution without taking the tidal force into account 



I force j 



The gas number density profiles of X-ray clusters of galaxies 
can be approximated by the empirical isothermal (3 model 
(Cavaliere & Fusco-Femiano 1976) 




n(r) 




-3/3/2 



(40) 



where no is the central gas number density, r c is the core 
radius of the cluster and (3 is the beta parameter. 

Recent high-resolution N-body simulations have sug- 
gested that dark halos of clusters are described by a family 
of universal density profiles. Navarro et al. (1997) proposed 
the following profile (NFW profile): 



0.2 



0.4 



. 6 



. 8 



Fig. 3. Functions F/{p\S) vs u (solid line) and G/(piS) vs 
u (dashed line). 



Because of the coincidence of these functions, the func- 
tion G(u) will be useful later. Thus, we can consider that 
the retarding function is proportional to the substructure 
velocity. 



Pdm{t) 



S c p cr QM 



(r/R s )(l + r/R s ) 2 
where the characteristic dimensionless density is 
200 c 3 



So 

R s 



3 ln(l + c) - c/(l + c) 



(41) 



(42) 



where c is a concentration parameter (Navarro et al. 1997). 
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T=t- 



It is useful to introduce a dimensionless time 



Rvirl 

The initial conditions x(0),v(0) are 
x(0) = -R vir i 
v{0) = -j-(0) = w 

(XT C s 



(43) 
(44) 



where v = dx/dr. 

The gravitational force acting on the X-ray gas of the 
subcluster equals 



F gr {R) 



F gr (R) 



GMg as z 1-2 , . , 

— ^2 — / ^nx p DM {x)dx 
4irGM gas2 5 c p cr £lM R 3 S 



R 



In 1 + — - 



R 2 

R/R s 



(45) 
(46) 



Rj l + R/R 6 



where M gaS 2 is the gas mass of the X-ray subcluster. 

The equation of motion of the X-ray gas of the substruc- 
ture without taking into account the tidal force (only the 
gravitational Eq. 46 and retarding forces Eq. 35 are taken 
into account) is then: 



d 2 x 20 n(x)kT S R v i r i dx 
dr 2 3 M gas2 c 2 dr 



4irG6 cPcr n M Rt x (47) 
i\ i i 



= 



The dependence of the X-ray substructure motion with 
time, which is obtained from Eq. (47) is shown in Fig. 
4 for the different values of the characteristic parameter 
(X/ujq = 0.25,0.5,1). 

The value of X/ujq = 0.5 corresponds to the first set of 
parameters, i.e., to 1E0657-56 (see Table 1). 




Fig. 4. Dependence of the X-ray substructure position 
x/R v on the dimensionless time tR v /c s for different char- 
acteristic parameters A/wo=0.25 (solid line), 0.5 (dashed 
line), 1.0 (dotted line). 

The characteristic parameter is useful for a classification 
of the X-ray substructure motion by analogy with an oscil- 
lator. If A/cjo < 1 the oscillator motion is a dying oscilla- 
tion, while if A/wo > 1 the oscillator motion is an aperiodic 
damping. 



The equation of motion of the dark (collisionless) matter 
substructure is (no retarding force): 



~j 2 -wGd c p cr \lMR s 2 3 X 

(It c s y 



(48) 



3.3. The characteristic parameter of the problem 

The second term on the left-hand side of Eq. (47) weakly 
depends on the space coordinate in the cluster core. When 
x -C R s , the third term on the left-hand side of Eq. (47) 
can be simplified. One can derive the Taylor series 



In 1 



R., 



Rs J Rs \ R 
The simplified expression of Eq. (47) is 



„ 20n o kTSR vir i , ATrGRf, irl S c p cr fl M n 
x H tttz — x H ^— s x = 



3M„ Qs2C 2 



2c? 



(50) 



This equation is the same as the harmonic oscillator one 

x" + 2Ax' + u 2 x = (51) 

The characteristic parameter of such an oscillator motion 
is 



A 



WmkTS 



3M gas2 c s V 2irG5 c p cr Vl 



1 



(52) 



3.4. Solution taking the tidal force into account 

It is necessary to consider the influence of the substruc- 
ture dark matter halo motion on the X-ray substructure 
motion. The assumption that the small cluster (the dark 
matter halo and the gas substructure) is a point mass is 
not acceptable because of the infinite gravitational energy 
at the initial time. It is a reasonable approximation to cal- 
culate the tidal force by assuming a NFW profile for the 
dark matter substructure halo. The tidal force acting on the 
subcluster dark halo is much smaller than the gravitational 
force between the two dark halos, because the tidal force is 
proportional to the mass of the X-ray gas of the subcluster. 
Consequently the dark matter substructure motion y(r) is 
described by Eq. (48). However, we must include the tidal 
force in Eq. (47) for the motion of the X-ray substructure 
by analogy with Eq. (46). 



Fudai — 47rM gas2 v £ 1 ScPctVLmR^-} 



In 1 



Rs 



Rs 



\ z ~y\ 

Rs2 



(53) 



where z is the position of the X-ray subcluster relative to 
the main cluster center in the case where the tidal force is 
taken into account. 

3.5. Mach number 

One task is to find the Mach number in the pre-shock re- 
gion. This is used to determine the temperature ratio of the 
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post-shock region to that in the pre-shock region. One can 
obtain the dependence of the Mach number in the pre-shock 
region on the Mach number of the piston: 



«2 



v sh -U (7 - l)Ml + 2 



Vsh 



(7 + mi 



(54) 



where v<i , v\ are velocities of the post-shock and pre-shock 
regions. Consequently, in the frame reference of the shock, 
v s h — c\M\ is the velocity of the shock, U is the velocity of 
the piston. 

With 7 = 5/3, one can find from Eq. (54): 



U 



3 Ml 
4 



1 



Mi 



(55) 



Therefore the Mach number of the piston is smaller than 
the Mach number in the pre-shock region. 



distance 



-0. 




Fig. 5. Dependence of the X-ray substructure position 
z/Rviri (dashed line) and of the dark halo position y/R V i r i 
(solid line) on the dimensionless time tR V i r \ / c s for the first 
set of parameters (1E0657-56). 



3.6. Examples: 1E0657-56 and Abell 1763 

The two sets of selected cluster parameters are in Table 1: 

1) For the first type of mergers when the substructure 
has passed through the main cluster core, the data on the 
galaxy cluster 1E0657-56 (the "bullet" cluster) are used 
(Markevitch, private communication) 

2) For the second type of mergers when the substructure 
has not gone beyond the main cluster center, as in the 
cluster of galaxies Abell 1763 (Lima Neto & Durret. 2007). 



Table 1. Sets of parameters for 1E0657-56 and Abell 
1763 





1E0657-56 


Abell 1763 


Cluster mass Mi, 10 15 Mq 


2 


0.5 


Central density no, (cm -3 ) 


0.005 


0.0079 


Core radius r c , (kpc) 


180 


136 


Mass ratio Mi/M 2 


60 


22 


Subcluster density, (cm~ a ) 


0.02 


0.0079 


Time unit Rvi r i/c s , 10 y j/r 


3.4 


2.8 




0.5 


3.1 



where the NFW concentration parameter is 4 in both 
cases. 



3.7. 1E0657-56 

The results can be compared with those of Clowe et al. 
(2006). The cluster exhibits a prominent bow shock with 
a Mach number M = 3.0 ± 0.4. The distance observed be- 
tween the two centers, corresponding to the X-ray and dark- 
matter components of the smaller subcluster, is around 200 
kpc, and the distance between the centers of the dark- 
matter substructures is around 1 Mpc. 

We approximate the shape of the X-ray subcluster with 
a cylinder (radius R=79 kpc and length L=120 kpc). For 
the first set of parameters, the position of the X-ray sub- 
structure is derived from generalized Eq. (47) taking the 
tidal force given in Eq. (53) into account, and the position 
of the dark halo is calculated from Eq. (48); these quantities 
are plotted in Fig. 5. 

The dimensionless time tR v /c s = 0.51 corresponds to 
the observed distance of 1 Mpc between the centers of the 
two dark-matter substructures and is used to determine the 



Mach number of the shock and the distance between the X- 
ray and dark-matter components. 

The Mach number of the X-ray substructure dz/(c s dt) 
is shown in Fig. 6 and the distance between the X-ray and 
dark-matter components y(r) — z(t) is shown in Fig. 7, as 
functions of the dimensionless time tR v /c s . 

From the condition tR v /c s = 0.51, we find the value of 
the Mach number of the X-ray substructure to be U/cx = 2 
from Fig. 6 and the distance between the X-ray and dark- 
matter components to be 200 kpc from Fig. 7. The Mach 
number of the shock Mi = 3 is found by solving the 
quadratic Eq. (55). 

Numerical hydrodynamical simulations of the bullet 
cluster 1E0657-558 (Springel & Farrar 2007; Milosavljevic 
et al. 2007) reached a very similar conclusion to ours: the 
halo collision velocity need not be the same as the inter- 
galactic gas shock velocity, and while the kinematics of the 
shock are sensitive to the details of the cluster structure, the 
instantaneous shock velocity can exceed the relative veloc- 
ity of CDM halos. 2D numerical simulations provide much 
more precise results than our approximate theory, because 
we made several simplifying assumptions (e.g., the problem 
is 1-D), and for accurate description of merging clusters 
numerical simulations are very useful. On the other hand, 
numerical simulations provide information for very specific 
sets of parameters. One immediate question is what are the 
key controlling parameters and might they affect the re- 
sults? For example, using different simulation parameters 
Milosavljevic et al. (2007) and Springel & Farrar (2007) 
find that the velocities of the CDM halos are respectively 
~16% and ^40% lower than that of the shock. Milosavljevic 
et al. (2007) strongly emphasize that the shock kinematics 
and morphology are extraordinarily sensitive to the param- 
eters of the simulations. The advantage of our approxima- 
tive theory is that it has a single characteristic parameter 
X/ujo, which completely determines the X-ray substructure 
motion. 

The similar results of Milosavljevic et al. (2007) and 
Springel & Farrar (2007), i.e., that the shock velocity is 
higher than the dark halo velocity, are expected for merging 
clusters, where the characteristic parameter X/loq is much 
smaller than 1. The influence of the retarding force in these 
systems is negligible and the velocity of the shock is higher 
than those of the X-ray substructure and the dark halo, 
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because the shock outdistances the piston in the "piston in 
a tube" model. 



Mach number 



3.25 
3 



1 . 75 




Fig. 6. Mach number of the X-ray structure vs the dimen- 
sionless time tR V i r \/c s for 1E0657-56. 



distance 




time 



Fig. 7. Dependence of the distance (kpc) between the X-ray 
and dark-matter substructures on the dimensionless time 
tRyiri/cs for 1E0657-56. 



3.8. Abell 1763 

We assume that the X-ray subcluster has the shape of a 
sphere and the radius of the X-ray subcluster equals the 
radius of the hottest region in the cluster, that is R w 
140 kpc. We propose that the hottest region in Abell 1763 
(where the temperature is twice the average temperature) 
is situated in the post-shock region. From the Rankinc- 
Hugoniot condition we can estimate the Mach number in 
the pre-shock region. The ratio of the temperatures in the 
two regions is: 

T 2 _ (2 7 M 1 2 -( 7 -l))(( 7 -l)M 1 2 + 2) 

T\ ( 7 + l) 2 M? (56) 

Using Eq. (55) we find the Mach number of the piston to 
be U/ci = 1.1. One can see from Fig. 9 that the dark 
matter substructure has a much higher velocity than the 
X-ray substructure. Therefore the effect of the retarding 
force is stronger in this cluster than in 1E0657-56. One 
can therefore expect a separation of the dark matter halo 
of the substructure and the X-ray substructure (Fig. 10). 
If the galaxies and X-ray subclump are not segregated, a 
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Fig. 8. Dependence of the X-ray substructure position 
z/Rviri (dashed line) and of the dark halo position y/ R V i r \ 
(solid line) on the dimensionless time tR V i r \jc s for the sec- 
ond set of parameters (Abell 1763). 

galaxy concentration near the region of the substructure 
dark matter halo is expected about 200 kpc (time = 0.6) 
from the center of the cluster Abell 1763 in the direction of 
the merger (south east to north west). 
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Fig. 9. Mach number of the piston (dotted line) and of the 
dark halo (solid line) vs dimensionless time tR vir \jc s for 
Abell 1763. 
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Fig. 10. Position of the dark matter substructure halo 
(kpc) vs dimensionless time tR V i r \/c s for Abell 1763. 

In order to check this point, we analyzed an optical im- 
age of Abell 1763 in the R band taken with the LFC camera 
on the Palomar telescope and kindly provided by D. Fadda 
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and A. Biviano. The SExtractor software was applied to 
the image, and stars were discarded on the basis of a plot 
of effective radius as a function of magnitude, leading to 
a catalogue of galaxies. From a magnitude histogram, the 
galaxy catalogue was estimated to be complete up to R 
~24. We therefore truncated the catalogue at this magni- 
tude to avoid effects due to incompleteness. 

Galaxies from this catalogue were then counted in 
squares of 140x140 kpc 2 centered on the cluster cen- 
ter. The numbers of galaxies in each square were found 
to be between 24 and 52, with the lowest number in 
the NW region, where our model predicts an overabun- 
dance of galaxies. The contamination by background galax- 
ies is estimated to be 9 galaxies in each square, as de- 
rived _frOTn_the_AAD_jnodel_l2^^1and-Hawthorn & Ellis 
( http : / /www. aao.gov.au/ astro/GalaxyCount/| . 

There is therefore no galaxy excess 200 kpc NW of the 
cluster center, indicating that the X-ray gas of the subclus- 
ter must have stopped while the galaxies have had time to 
travel further out. Hence, a direct correspondence between 
the hottest region in Abell 1763 and the post-shock region 
is not correct. Abell 1763 is therefore a structure in which 
the observed baryons of the subcluster and the galaxy con- 
centration are spatially segregated; a full analysis of the op- 
tical image will be done elsewhere and is beyond the scope 
of this paper. 

Dark matter annihilation in the halos of galaxy clusters 
has relevant astrophysical implications (e.g., Colafrancesco 
2006). In fact, if dark matter is constituted by weakly in- 
teracting massive particles (WIMPs), their annihilation can 
produce direct observable signals. However, the spatial and 
spectral intensities of the astrophysical signals coming from 
XX annihilation are expected to be confused with, or even 
overcome by, other astrophysical signals originating from 
the intracluster gas and/or from the relativistic plasmas 
present in the cluster atmospheres, especially when all these 
components are co-spatially distributed with the DM com- 
ponent. An ideal system to detect a DM annihilation signal 
would therefore be a cluster with a clear spatial separation 
between the various matter components (Colafrancesco et 
al. 2007). 

4. Evaluation of the incidence of substructures 

Clusters of galaxies form by the mergers of smaller clus- 
ters and groups. The presence of substructure in clusters of 
galaxies can be used to estimate the cosmological density 
parameter Qm- 

X-ray data allow the quantification of substructures be- 
cause cluster mergers compress and heat the intracluster 
gas, and this can be measured as distortions of the spatial 
distribution of X-ray surface brightness and temperature 
(Forman & Jones 1994; Jones & Forman 1999). An alterna- 
tive method at optical wavelengths to study substructures 
is based on the distribution of cluster galaxies. However, the 
optical method requires a large number of galaxies (Dutta 
1995), at least a few hundred, and is more susceptible to 
contamination from foreground and background objects. 

X-ray studies of cluster substructure use a number 
of different statistics (see Buote 2002 for a review). 
Several studies have used the ellipticity (Gomez et al. 
1997; Kolokotronis et al. 2001; Melott et al. 2001; Plionis 
2002), while other studies measured centroid (or center- 
of-mass) variations (Mohr et al. 1995; Gomez et al. 1997; 



Kolokotronis et al. 2001). Another method developed by 
Buote & Tsai (1995) is the power ratio method, which is 
constructed from the moments of the X-ray surface bright- 
ness. 

Richstone et al. (1992) performed the first theoreti- 
cal study of the relationship of substructure to cosmology. 
They assumed that any substructure is wiped out in a clus- 
ter crossing time and calculated the fraction of clusters in 
the spherical growth approximation that formed within the 
last crossing time as a function of Q m and A. This frac- 
tion primarily depends on J7m, and they estimated that 
£Im > 0.5 based on estimates of the frequency of substruc- 
ture in low-redshift clusters from X-ray images. 

The approach proposed by Richstone et al. (1992) to 
determine VLm suffers several limitations. In their semi- 
analytical calculations they avoided the issue of the power 
spectrum by concentrating on clusters having the same to- 
tal mass. The collapse time of a 1O 15 M0 spherical den- 
sity perturbation was defined to be the dividing point be- 
tween clusters that do and do not possess substructure. The 
relationship between the collapse time of a spherical den- 
sity perturbation and subclustering, though qualitatively 
reasonable, is ambiguous. Consequently, they predict an 
ambiguous "frequency of substructure" rather than a well- 
defined quantitative measure of cluster morphology. 

A more recent semi-analytical approach was developed 
by Buote (1998). He assumed that the amount of substruc- 
ture depends on the amount of mass accreted by a clus- 
ter over a relaxation timescale and relates the fractional 
accreted mass to the power ratios. Although these semi- 
analytical methods give an indication of the expected evo- 
lution of cluster substructure and its dependence on cosmo- 
logical parameters, perhaps the best method to constrain 
cosmological models is through the comparison to cluster 
simulations. 

Evrard et al. (1993) employed N-body/3D gas dynamic 
simulations of the formation of galaxy clusters and demon- 
strated the existence of a morphology-cosmology connec- 
tion for X-ray clusters. Confirming the analytical expecta- 
tions of Richstone et al. (1992), they showed that cluster 
evolution is sensitive to the cosmological model in which 
the clusters form. Numerical simulations show that both 
the centroid shift and the power ratios are capable of dis- 
tinguishing cosmological models (Evrard et al. 1993; Buote 
& Xu 1997; Valdarnini et al. 1999; Suwa et al. 2003). Buote 
& Xu (1997) find that the power ratios of their ROSAT 
clusters indicate an Qm < 1 universe. Valdarnini et al. 
(1999) find that the ACDM model is inconsistent with the 
data. However, they used as = 1-1, which is fairly high and 
may cause the disagreement between the ACDM model and 
the data. Suwa et al. (2003) find that the center shifts and 
power ratios of both surface brightness and projected mass 
density are able to discriminate between the cosmological 
models at z=0. 

We uses the value of Qm — 0.268 ± 0.018, which was 
obtained from the WMAP three year observations (Spergel 
et al. 2007), to find the dynamical time of the substruc- 
ture, which is higher than the gravitational free-fall time. 
A plausible reason for this disagreement is that the fraction 
of clusters with X-ray substructures exceed the number of 
clusters that have collapsed during the last gravitational 
free-fall time due to the retarding force acting on the X-ray 
substructure. 
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4.1. The collapse time for bound perturbations 

The separation of two observers in a Friedman-Robertson- 
Walker universe with a cosmology constant obeys the equa- 
tion of motion 

d 2 r AnGp+ij Ar 

^2= ^ + T (57) 

where , r* are the mean density of the universe and r is 
the separation between two observers at any specified epoch 
in the universe. 

This equation also describes the evolution of the central 
distance of a point anywhere in a homogeneous perturba- 
tion. 

The equation of mass conservation shows that £Im{z) 
and A(z) of the universe evolve as 

(i + z) 3 n M0 



X(z) 



H 2 {z)/m 

Ao 



(58) 



H 2 (z)/H$ 

where fi = 8nGp/{3H 2 ), A = A/{3H 2 ). 

The dependence of the rate of expansion H(z) is given 

by 

- (l-r!A./o-Ao)(l + ^) 2 + r! M o(l + 2) 3 + Ao (59) 



One can calculate from Eq. (65) the critical value of per- 
turbations 6 cr = 2.5 • 10 -3 . 

The maximum expansion factor for a homogeneous spheri- 
cal perturbation is given by Eq. (62). 

Direct substitution in Eqs. (60), (61) and (63) yields the 
ratio of the collapse time to the current age of the (unper- 
turbed) universe for bound perturbations 

t _ 2 f[ m/ri du (i -n M {i + 5) - x + n M (i + S)^ 1 + \u 2 )~ 1/2 
To ~ j 1+zi du (i - n M - x + Qmu- 1 + Xu 2 y 1/2 

where £Im, X are fixed at the values evaluated at 1 + z\. 

4.3. The estimation of the dynamical time 

For a particular universe, specified by S1mo,Ao, Eq. (66) 
defines a monotonic relation between S and t/Tq for bound 
perturbations, and can therefore be regarded as an equation 
for S as a function of t/Tq (Richstone et al. 1992). Eqs. (64) 
and (66) specify, for a particular fij^O) Ao, cr, the fraction of 
the universe on some specific mass scale that has already 
collapsed at cosmic time t: 



Hi 



F(t/To) 



It is possible to rewrite Eq. (57) in integral form and esti- 
mate the time to expand from r+ to r for a homogeneous 
perturbation: 



T = H' 1 



r fr* 



fi - n M * - x* + n+u 1 + x*u 2 ) 1/2 



du(60) 



Eq. (60) is only valid for perturbations that expand at least 
as far as r/r*. Bound perturbations which do not expand 
forever reach a maximum expansion and recollapse on a 
time scale given by 

r = 2T{r m /n) (61) 
where r m is the first real root > 1 of the equation: 
(1 - n Mi - A*) + flu*™- 1 + X*u 2 = (62) 



4.2. Initial density perturbations 

Now suppose that at the epoch defined by 1 + Z\ = 1000 
the universe is characterized by an unperturbed Hubble 
flow and cosmological constant, but that the density, on the 
mass scales of clusters of galaxies, is normally distributed 
about the mean density with variance a 2 . Thus, we can 
write for the perturbations 5: 



A* = X( Zl ) 

^m* = ^m(zi) ■ (1 + 5) 
and 5 has the Gaussian distribution 



dF(S) = 



1 



V2 



: exp 



7T(J 



2a< 



dd 



(63) 



(64) 



The bound perturbation of the density must exceed the 
critical value S cr , which is given by 

l-n M (zi)-(l + 5 cr )-X(zi)+ (65) 

+1.5 ■ [2A(zi) ■ (fijif ■ (1 + Scr)) 2 } 1 / 3 - 



1 - erf 



( S(t/T ) 
V V2a 



(67) 



where erf is the error function. 

Eq. (67) predicts the fraction of the universe currently 
in virialized clusters with a mass scale of rich Abell clusters. 
That fraction is the product of the number density < n > 
and mass of rich clusters M, divided by the mean density 
of the universe. The value of a must be found by solving 
Eq. (68): 



1 -(l-er f m)) = <n>M 

2 V \V2aJJ Pcr^MO 



(68) 



The simplification of Eq. (66) is useful from an analytical 
point of view in order to find the function 6(t/Tq). A good 
approximation can be found by assuming 



-(5) oc 8- 3 ' 2 



(69) 



The function t/Tq computed from Eq. (66) is shown in Fig. 
11. 
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Fig. 11. Collapse time t/Tq as a function of the density 
perturbation 5. 
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The fraction of the mass of present day clusters that 
have collapsed within the time interval (To — 5t, To) given 
in Fig. 12 is 



(70) 




interval 



Fig. 12. The fraction of the mass of present day clusters 
that have collapsed within the last time interval St/To. 

Forman & Jones (1994) have studied the X-ray sur- 
face brightness of 250 clusters of galaxies observed with 
Einstein. They concluded that 30% of that sample had X- 
ray substructure. Taking into account this observational 
constraint one can derive the dynamical time, and we find 
St /T = 0.17 (the big point in Fig. 12). 

4.4. The retarding of an X-ray substructure 

Since Eq. (51) can be solved analytically, we estimate the 
retarding of a X-ray substructure from this equation when 
the characteristic parameter A 2 /uj 2 >> 1. 
The solution of Eq. (51) in case X/ujq > 1 with initial con- 
ditions (x(0) = -R,v(0) = 0) is 



x = (Ci + C 2 ) exp(-(A - VA 2 -uj 2 ) ■ t) + 
+ (d - C 2 ) exp(-(A + VA 2 -^ 2 ) • t) 
where 



d = -R/2 
<h = - 



(71) 
(72) 

(73) 



XR 



2^/X 2 



UJZ 



The X-ray substructure velocity is found from the Eq. (71): 
dx 



^ = Ae- Xt ch{^\ 2 - wg • t) + Be- xt sh{^JX 2 - uj 2 • t)(74) 
where 



A = 2{-Xd + J X 2 -lo 2 C 2 ) 



(75) 



B = 2{-XC 2 + sJX 2 -u 2 C 1 ) 
If the characteristic parameter X 2 /u}q 3> 1, we have 

,2 



exp(-At + \/X 2 - ojfit) w exp ( ~^r* 



(76) 



Therefore when the time t 3> 1/A, we can consider only 
the terms of Eq. (73) that contain the factor exp(— At + 
y/X 2 — uJqI). The dependence of the velocity with time (t >• 
1/A) is 



— w i?^exp --£f 
dt 2A V 2A 



(77) 



exp(-At - VA 2 - ujH) w exp(-2Ai) 



The derived velocity of the X-ray substructure in case 
A/ wo > 1 is much smaller than the dark halo velocity 
(~ Ruiq). Therefore for merger events, when this character- 
istic parameter is much higher than 1, the dynamical time 
is much higher than the gravitational free-fall time. This ar- 
gument can explain that the merger duration is higher than 
the gravitational free-fall time for many merger events. 

5. Conclusions 

There are many proofs of evidence that X-ray substructures 
move with a supersonic velocity in the ICM of clusters. 
In their recent review of shocks and cold fronts in galaxy 
clusters, Markevitch & Vikhlinin (2007) nicely summarize 
our present knowledge on physical mechanisms taking place 
when clusters merge, such as those discussed in the present 
paper. The analytical study of the X-ray substructure mo- 
tion is proposed in this article. The retarding force must 
be taken into account to obtain a correct solution. One of 
the most obvious manifestations of its importance is the 
spatial segregation of the stellar component and the X-ray 
substructure in the galaxy cluster 1E0657-558. The behav- 
ior of the X-ray substructure depends on the characteristic 
parameter (see Fig. 4) in the "piston in a tube" model. 
When the characteristic parameter is smaller than 1, the 
X-ray substructure can pass through the main cluster core 
(e.g. 1E0657-558). We also examine the case when the char- 
acteristic parameter is higher than 1 and the effect of the 
retarding force is crucial (e.g. Abell 1763). The shock ve- 
locity in front of the X-ray substructure is calculated. The 
method for predicting the position of the dark matter halo 
in clusters, where the observed baryons and the inferred 
dark matter are spatially segregated, is considered. 

Richstone et al. (1992) used the dynamical time of 
the substructure to constrain Qmj but their constraint 
fi&f > 0.5 does not correspond to the current data. A pos- 
sible reason for this is that substructure lasts significantly 
longer than a crossing-time in clusters. The retarding force 
influences the substructure dynamical time, which is then 
higher than the crossing-time. Consequently, an accurate 
description of substructure propagation is a requisite to es- 
timate SI a/. 
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